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The quantification and classification of quantum entanglement is a very important 
and still open question of quantum information theory. In this paper, we describe 
. an entanglement measure for multipartite pure states (the minimum of Shannon's 

I entropy of orthogonal measurements). This measure is additive, monotone under 

LOCC, and coincides with the reduced von Neumann entropy on bipartite states. 
' A method for numerical calculation of this measure by genetic algorithms is also 

presented. Moreover, the minimization of entropy technique is extended to fermionic 
states. 
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=3 : I. INTRODUCTION 

Multipartite quantum entanglement plays significant role in whole quantum science. For 
example, one can easily obtain that we can't achieve algorithmic speedup in quantum com- 
puter without quantum entanglement. Also entanglement is needed for other quantum 
I protocols: quantum teleportation, quantum error correction, etc. What is more, the theory 

CN ■ of quantum entanglement may help us to understand multiparticle quantum physics deeply. 

. It's well known that the entanglement of bipartite pure states is fully described. Schmidt 

^ i coefficients uniquely define a local unitary orbit of a bipartite quantum state and unam- 
Q>^ I biguously determine the class of states that can be obtained from by LOCC [if. 
O ■ Schmidt decomposition 



® \j)B = v^K)a ® \i)B 



is an equivalent of matrix SVD decomposition 

E Cij\i)A{j\B = E y/>H\i)A{i\B, 



A = USV*. 



It's easy to see that there is no Schmidt decomposition even for three qubits. The W 
state 

-^(|001) + |010) + |100)) 

cannot be represented as 

Ao|0) ® |0) ® |0) + Ai|T) ® |T) ® |T). 



* Electronic address: | andrey. chernyavskiy @gmail .com 
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A stronger counterexample exists for SVD decomposition: there is no lower-rank orthogonal 
approximation even for third-order tensors 

While a lot of applications use SVD: image compression (for example, LSI (Latent 
Semantic Indexing, for example, [3]), statistics (for example, jsj), etc., the appropriate 
expanding of SVD decomposition (Schmidt decomposition) to higher order tensors (to more 
than two subsystems) is an important task not only for quantum information theory but for 
fundamental (tensor algebra) and applied (statistics, machine learning, etc.) math. Some 
higher order analogues of SVD may be found (for example, |6|), but non of them can be 
directly applied to the quantum pure state case. 

Whereas there is no full theory of multipartite quantum entanglement, some important 
results are known. 

A full classification of pure states of three qubits in terms of SLOCC (stochastic LOCC) 
is given in [7|. A description of quantum entanglement with nilpotent polynomials is con- 




Some entanglement measures of pure and mixed states can be found, for 
Theory of quantum entanglement measures is studied in. 



16, 17, 18, 19 



structed in 
example, in 
for example 

The paper is organized as follows: In the first section of this paper we present the min- 
imum of Shannon's entropy of orthogonal measurements and prove that this function pos- 
sesses all necessary pure state entanglement measure properties. Some other features of 
this measure and its values for some states are also presented. Next we describe a method 
to compute this measure using genetic algorithms. An extension of minimal measurements 
entropy to fermionic states is described in the last part of the paper. 



A. Notation 



Hsh is the Shannon entropy, 
H^N is the von Neumann entropy. 



II. SHANNON'S ENTROPY OF ORTHOGONAL MEASUREMENTS 

MINIMUM 

Consider a pure n-qudit qauntum state 

d 

IV') = «n,i2,...,i„Kl^2 •••««)• (1) 

il, 12,..., in = l 

The measurement of this state in the computational basis gives states \i1i2 ■ ■ - in) with 
probabilities |aij,i2,...,i„p. So we can regard this measurement process as a signal generator. 
We shall say that Shannon's entropy of this generator is called measurements entropy of the 
state \ip): 

d 

Hmeas{\ij)) = Hsh{Diag{\i)){'4)\)) = ^ \ai^,i2,...,iS- 

jl,«2,...,in = l 

But this characteristic of a quantum state is not invariant under local changes of the mea- 
surement basis. Shannon's source coding theorem shows that, in the limit, the average 
length of the shortest possible representation to encode the messages in a given alphabet 
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is their entropy divided by the logarithm of the number of symbols in the target alphabet. 
That's why the natural way to construct an invariant is the minimization of -ffmeasd"^)) over 
all possible local changes of the measurements basis: 



and Enmin is a measure of entanglement of the pure state \ip)- 

Remark II. 1. The definition of EHmm for subsystems with different dimensions is the same 
to 1^. As non of the following reasoning hears on the equality of subsystem dimensions, here 
and further the common d for different ij in (J\) is used only for notation simplicity. 



Entanglement measure for pure quantum states must satisfy the following intuitive con- 
ditions: 

(i) must be equal to zero for fully unentangled states 

(ii) must be invariant under local unitary operations 

(iii) must be invariant under attachment and detachment of unentangled ancilla 

(iv) must not increase under LOCC. 

Let's show that Eumin satisfies (i)-(iv). First of all we consider another one important 
property: 

Property II. 2. (Additivity of EHmm) 

Let and \il)2) be multipartite (may he entangled) pure states, then 



The proof is straightforward from the definition of EHmin and additivity of Shannon's 
entropy. 

The properties (i)-(ii) are trivial and (iii) is the consequence of EHmin additivity. To prove 
(iv), since (ii) and (iii) are correct, we may only prove that EHmin doesn't increase under or- 
thogonal measurements in the mean. To show monotonicity under orthogonal measurements 
in the mean we need some lemmas. 

Lemma II. 3. Let pab be a density matrix of a pure bipartite state. Then 

H,h{Diag{pA)) < H,h{Diag{pAB)) < Hsh{Diag{pA)) + Hsh{Diag{pB)), 
where pA and ps are reduced density matrices. 

Proof. The left inequality follows from generalized grouping of Shannon's entropy: 



Enmin 



min Hmeas{Ui ® U2 IS) ■ ■ ■ IS) Unlip)); 



(2) 



A. Necessary properties of entanglement measure 



Enminili^l) ® IV'2)) = EHminili^l)) + E Hmini\'4^2)) ■ 



HshiPl, ■■■ 



Pct„_i + 1) ■••) Pan 




i=l 



j=o-i_i + l 



j=(7i-l + l 
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The right part follows the sub-additivity of entropy: 

n n n m m m 

i=l i=l i=l j=l j=l j=l 



Lemma II. 4. Let 

d 



□ 



1^) = «n,i2,...,inl^l^2 

be a pure n-qudit state. The states \j) ® \ipj) are results of the measurement of in 
computational basis with probabilities pj. Then 

d 

Hmeas{Ul\lp)) > Pj Hmeas{\j) ® IV'j)), (3) 

where Ui is an arbitrary unitary transformation of the first qudit. 

Proof. If |j) ® \ipj) are results of measurement in computational basis with probabilities 
Pj, then can be represented in the view 

d 



where 



\cj\^ =Pr 



Let ipj,k = 1, be the amplitudes of l^/'j), then 

d 

= H,,{\cim\^ • • • , icii vr?, ■ ■ ■ , iqi v^r, • • • , ic.ri^r?) 

= { by the strong additivity of Shannon's entropy } 



= J2pjHsh{\^P]\\ |^,T, . . . , T) + H,h{pi,P2, ...,Pd) 

i=i 

d 

Using this we can rewrite in the equivalent form 

Hmeasim < Hmeas{Uim + H,h{p). (4) 
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Let p = be the density matrix of the state Pi = r'"2,3,...,d(p), P2 = Tri{p). 

Then by Lemma IIL3I for p we get 

Hsh{Diag{p)) < H sh{Diag{p2)) + Hsh{Diag{pi)), 

or equivalently 

H^eas{m < Hsh{Diag{p2)) + Hsh{v)- (5) 

Let pui = Ui\il))l{il)\Ul be the density matrix of the state Ui\il)), p^^ = Tri{pu^). Since 
Ui affects the first qudit only, p^^ = p2- 
By Lemma [11.31 

H,h{Diag{p2)) = H sh{Diag{p2^)) < Hsh{Diag{pu^)). (6) 

By ([5]) and ([6]) it follows that (j4]) is correct, consequently ([3]) is correct too. □ 

Theorem II. 5. (the monotonicity of EHmin under orthogonal measurements) 

Let lip) be an n-qudit state. The states are results of some orthogonal measurement 
with probabilities pj . Then 

T^Hmin 

j 

Proof. We may assume, without loss of generality, that we measure first qudit in the basis 
{ipl} (a measurement of more than one qudit can be replaced by sequential one-qudit mea- 
surements). Then {ipj) = {ipj) l"^])- Let \ipmin) has minimal measurements entropy over 
local unitary orbit of \ip), and \ipmin) = Ui ® U2 ® ■ ■ ■ ® Unl^j) ■ Then 

^PjEHmin{\ijj)) = ^PjEHminili^]) ® < ^PjHmeas{\j) ® (f/2 ® • • • ® Un)\^pj)) 

j j j 

< {by Lemma Ullf < Hmeasii^min) = Enminili'))- □ 

Thus, from III. 5] we get (iv). 

B. Other properties 

Lemma II. 6. (Klein's lemma) 
Let p be a density matrix, then 

Hsh{Diag{p)) > H^n{p)- 

Theorem II. 7. EHmin coincides with the reduced von Neumann entropy for bipartite states. 
I.e., let \ip) be a pure bipartite state; p = is its density matrix, pa = Ttb^p) and 

Pb = Tta^p) are reduced density matrices of subsystems. Then EHrnin{\4')) = Hyj^{pA) = 
Hvn{pb)- 
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Proof. Consider the Schmidt decomposition of lip): 

IV") = a/A^Ka) ® lis)- 

i 

Then 

EHmin{\lp)) = min Hmeas{UA ® UbHj)) < HshiK) = H^n{Pa) = H^n{pb). 

From Lemma III.3I we get 

Hmeasil'ip)) = Hsh{Diag{p)) > Hsh{Diag{pA))- 

Further, by Klein's lemma we have 

Hsh{Diag{pA)) > H^n{pa) = H^n{pb)- 

Thus EHmin{\'4')) = H,,n{pa) = H^n{pb)- □ 
Property II. 8. (EHmin of generalized GHZ states) Consider a generalized GHZ-state: 



\GHZ) = ^ai\i)-^ ® \i)2 ® K)„, 



i=l 

d 

where G M, ^ |ajp = 1. 

i=l 

Then 

Hmeas{\GHZ)) = miu Hmeas{Ul ^ U2 ^ . . . ^ Un\GH Z)) . 

Ul,U2,...,Un 

Proof. We can consider the space of GHZ state as bipartite space of the first qudit and the 
remaining qudits. Then 

d 

\GHZ) = ^ai\i)^ ® 

i=l 

From Theorem III.7I and Remark III.ll we have that Schmidt decomposition has minimal 
measurements entropy. From this we get 

d 

HmeasilGHZ)) = H^easiJ^ ^^\') 1 ® K)) ^ Hmeas{Ul ® U\GHZ)), 

i=l 

where Ui is a unitary transformation of the first qudit and U is an arbitrary unitary trans- 
formation of 2, 3, ■ ■ ■ ,n qudits. We can take ® ® ■ ■ • C?) f/„ as f/, then 

this finishes the proof. □ 

d 

So we have EHmin{\GHZ)) = - XI |ajpln |aip. 

2=1 
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C. Numerical properties 

The following properties were obtained numerically by genetic algorithms. 

Numerical Result II. 9. (EHmin of generalized W states) 
Consider a generalized W state 

\W) = ai|0 . . . 01) + Oslo . . . 10) + a„|l . . . 00), 

n 

where |ai| =1. Then 

i=l 

n 

EnminiW) = Hmeas{W) = ^ |aj|^ln \aif . 

1 

Numerical Result 11.10. Let 

d 

n,j2,.--,2n = l 

d 

I'f) = K,i2,...,in\hi2---in), 

il,i2,...,in=l 

and Hmeasili')) = HmeasiW)) = EHmin{\i^)) = EHmin{W)) (^^ otker words and \ip) have 
equal EHmin o.iT'd they are in minimal entropy representation), then 

la- ■ ■ |2 = 16. . . |2 

to within local permutations of basis vectors. 

I. e. modulus squares of a minimal entropy representation of the local unitary orbit are 
unique. 

D- EHmin IS Substantially multipartite 

Numerical computations using genetic algorithms show that two equivalent under bipar- 
tite entanglement states of three qubits \ip) and I?/') may have different EHmin- This means 
that EHmin IS Substantially multipartite. Equivalence under bipartite entanglement means 

3r/^ r/^ jp IP IP IP ■ 

where f/f are unitary evolutions of i-th qubit, and U^j are unitary evolutions (may be 
entangled) of i-th and j-th qubit. 

Prom EHminili^)) 7^ Enminil'p)) have that \(p) and are not equivalent under local 
unitary transformations, i.e 



$Ui, U2, Us : \ip) = C/i C/2 C/slV')- 
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This property of EHmin gives an answer to the problem of equivalence of bipartite and 
multipartite entanglement proposed in [20^1: the bipartite and multipartite entanglements 
are not equivalent. This also means that entanglement measures based only on Schmidt 
coefficients of different decompositions of a state are not good for quantifying multipartite 
entanglement. 



III. CALCULATION OF Enmin USING GENETIC ALGORITHMS 
A. Formalization of the optimization problem 

Consider the n-qudit state 



11,12, .■■,«« = ! 

Then 

{\^)) = min Hmeas{Ul ® f/2 (S) . . . ® Un\ip)), 
Ui,U2,-.-,U„ 

where f/j is a ci-dimension unitary operator on i-th qudit, 

d 



'll,l2,...,t-n 



To solve this optimization problem we need to parameterize unitary matrices f/j. When 
d = 2 this parametrization is well known: 

. , /e^(-'3'-''')cos7i -e^(-^»+*')sin7,- \ , . 

^^^^^'^^'^^^ - ^ e^(ft-'^»)sin7i e^(ft+'^>)cos7, j ' 

where Pi,6i, 'ji are real numbers. 

Parametrization of unitary matrices for d > 2 is a still open and interesting question. The 
best known parametrization is proposed in 21, 23], but it is very slow, so we have chosen 



another one. Let's parameterize hermitian matrix H by real numbers (the diagonal has 
d real numbers, and d{d — l)/2 complex numbers above the diagonal need d{d — 1) real 
numbers to be parameterized). Then we take U = e^^ as a unitary matrix. 

In such a way the calculation of EHmin is an optimization task of Sra real parameters 
for qubits and nd'^ real parameters for qudits {d > 2). To calculate -E//mm(|^)) we need to 
minimize the function 



r{Xi, . . ..Xn-k) = Hmeas{U{Xi, . . . , Xk) ^ U {x k+1, ■ ■ .,X2k)®-- ■ ®U {X(n~l)k+1, ■ ■ • , a;n-fc) | V^) ) , 

where U{xi, . . . , Xk) is a parametrization of d x d unitary matrix by k real parameters. 

In the general case, functions f^{xi, . . . , Xn-k) may be multimodal (have many local mini- 
mums), because of this fact we can't use gradient-based optimization methods. Thus genetic 
algorithm has been chosen. 
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B. Genetic algorithm (GA) 

There are a lot of publications about GA, but for a brief overview and references we 
recommend Wikipedia [23]. 

GA has been already used for quantum entanglement calculation (to calculate relative 
entropy of entanglement of mixed bipartite states) in [24 . 

Now we describe GA that was used for Enmin calculations. 

Every parameter x of is being encoded by n^en real value genes gj,j = l,ngen by the 

rule X = 10^ ■'fi'j- So, if we have k parameters, then a chromosome is a vector {gi} of 
i=i 

length Ugen ■ k. (A full tuple of parameters of is encoded by a chromosome.) 
A mutation of a chromosome {gi} is determined by the following probabilities: 

P{9T = 9i) = (1 - Pmut), P{gT = 9i + C)= Pmut, 

where {g^^} is a chromosome after mutation, Pmut is a mutation probability, ^ is a random 
variable that has uniform distribution on [—rrimuty'^mut]- 

Crossover of chromosomes {g}} and {gf} is a chromosome {g^}, where 

A fitness function is —f^- 
The algorithm. 

1. Initialize the first population of UpopuiaUon random chromosomes with uniformly dis- 
tributed on [-minit,minit] genes. 

2. Wait an epoch (this step will be described further). After the epoch we have a new 
population. 

3. If one of termination conditions is satisfied, algorithm stops and the fitness of the best 
chromosome of the last population is taken as a result, else we repeat from Step 2. 



Termination conditions: 

We can fix the maximal number of epochs(r;,epoc/is), the precision e, and the maximal 
number of nonchanging epoches Uterm- Then termination conditions will be: 
-We reach n^pochs epoche. 

-Examine the best chromosome from each of the last rtterm epoches. If the fitness function 
values of these chromosomes differ from each other by less than e, then this termination 
condition is true. 
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The epoch. 

1. Two random pairs of chromosomes are chosen from the best Upopuiation — nbad chro- 
mosomes of the population, Ubad is a number of the weakest chromosomes of the 
population, which cannot be used for reproduction. From each pair we take a chro- 
mosome with the best fitness. After that we put a crossover result of the two selected 
chromosomes to the new population. 

2. Repeat Step 1 while the new population is lesser than npopuiation. 

3. Mutate every chromosome of the new population. 

While for lip) from local orbits of "easy" states like GHZ or W has no "difficult" local 
minimums, the situation is reverse for states with random amplitudes. For example, one 
7-qubit state has a local minimum 4.0220, while another one (probably global) minimum is 
3.968. To solve this problem we use GA with "islands": we form Uisiands islands with rather 
small equal populations. Epoches at these islands are independent and only infrequent 
migrations are allowed. The increase of islands count is an equivalent of the the whole GA 
repeating, so the probability of an error decreases exponentially with the number of islands. 
Increase of the probability of migration speedups the convergence, but it also increases the 
probability of local minimums results. 

Remark (about software implementation) 

One of the advantages of GA is its parallelism, so the software realization of Enmin 
calculation was multithreading. Moreover, the implementation of EHmin calculation for 
qubits was realized using nVidia CUDA [25| technology. Using cheap personal GPU, 7x 
speedup against the best 4-core CPU implementation has been achieved. Software complex 
can easily calculate Eurnin for up to 17 qubits inclusively. 



IV. EHmin FOR FERMIONIC STATES 

Consider a pure state of n fermions in p-dimensional Hilbert space with basis /i, /2, . . . , /„: 



I/) 



E 

il,i2v;in = l,n<i2<---<«r7 



\ili2 ■ ■ . in), 



where 



lA |2 - 1 

il,j2,.--,«n = l,n<«2<---<iji 

fhin) ... /ii(r„) 



\i1i2 . . .in) 



1 



are Slater determinants. 

This states correspond to normalized elements of the exterior algebra A'^C'^. Let's take 
states that correspond to separable elements xi /\x2 /\. . . /\Xn oi A'^C'^ as unentangled, i.e. 
the states that are represented as a single Slater determinant in some basis of one-particle 
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space are unentangled. Thus, we can take unitary changes of basis as unentanghng 
transforms. 

This is more mathematical than physical approach. For example, if we have two electrons 
in two quantum dots their one-particle space must at the least include spin and position 
coordinates. In the simplest case the whole space will be H spin® H position and its basis will be 
I I) ® I I) ® |2), I I) ® I I) ® |2). As we can see, unitary transformations of this basis 
are "entangled" in an intuitive physical way. So, the construction of a physical formalism of 
unentangled local transformations of identical particles is an important direction for future 
research. 

But mathematical approach is nevertheless very important and is used in many papers 
devoted to the entanglement of indistinguishable particles. For example, the Slater decom- 



position (an analogue of the Schmidt decomposition) of two fermions is constructed in [26 



and we will discuss it further. Another interesting result is a classification of A^C^ fermionic 



states in terms of SLOCC 27 



Now let's define Efjmin for fermionic states. 

The change of basis is determined by the unitary matrix U and in the new basis the state 
becomes 



Uo\f)= J2 K^^2...^J^1^2■■.^n), 

ii,i2,...,in=i,ii<i2<---<in 



il,i2,.--,*n=l,il<i2<---<ip 

where Mf^'/^'"'f" is a determinant of the matrix that is constructed by the crossing of 
3i-,j2i ■ ■ ■ lin columns and ii,i2, ■ ■ ■ ,in rows of U. 
Measurements entropy of |/) is 

p 

Hmeasilf)) = ^ ^ | ■^jij2---«n I ^^g | Aj j j2...jjj | . 

il,j2,.--i*n = l,H<«2<---<in 

And 

EHminilf)) = rain HmeasiU O |/)). 

The parametrization of U is described in Section [III A[ To get U o \ f) amplitudes we 
need to calculate all minors of U, this can be done recursively using determinant expansion 
by minors. 



A. Slater decomposition 

Now consider a state of 2 fermions in 2p-dimensional space: 

2p 

1/2)= ^ Ki2\Hi2)- 

«1,«2 = 1,«1<«2 

Then this state can be represented as a linear combination of only p Slater determinants 
using a unitary change of basis [26| : 

p 

Uslater O I/2) = ^ Zi\2i, 2i + I) . 
i=l 
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This representation is called Slater decomposition. 

This fact follows from the existence of the symplectic basis for antisymmetric matrices 



28( 1 ■ full proof may be found in [26j . 

An important numerical result has been achieved for Slater decomposition and EHmin 

Numerical Result IV. 1. Let I/2) be the state of 2 fermions with 2p- dimensional one- 
particle space. Then the following two conditions are equivalent: 
1. The state I/2) is in the Slater decomposition representation: 



|/2) = 5^;2.|2^,2z + l). 



2- EHmin{\f2)) — -f^meas ( | /2 ) ) ■ 

(I.e. Slater decomposition, like Schmidt decomposition, has minimal measurements en- 
tropy, and the process of finding EHmin for a state of two fermions is equal to finding its 
Slater decomposition.) 



V. CONCLUSION 

Entanglement measure EHmin for multipartite pure states was presented. Also we proved 
that this measure is additive, satisfies all necessary entanglement measure conditions, and 
coincides with the reduced von Neumann entropy for bipartite states. The method of numer- 
ical calculation of EHmin by genetic algorithm was presented and tested on up to 17 qubits 
inclusively. Moreover, Enmin was generalized to fermionic states, and this generalization 
corresponds to Slater decomposition for two-fermions states. 
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